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We consider a classically chaotic system that is described 
by a Hamiltonian H(Q,P;x), where (Q,P) describes a par- 
ticle moving inside a cavity, and x controls a deformation of 
the boundary. The quantum-eigenstates of the system are 
\n(x)). We describe how the parametric kernel P(n\m) = 
|(n(a;)|m(:ro))| 2 , also known as the local density of states, 
evolves as a function of 8x = x— xq. We illuminate the non- 
unitary nature of this parametric evolution, the emergence of 
non-perturbative features, the final non-universal saturation, 
and the limitations of random- wave considerations. The para- 
metric evolution is demonstrated numerically for two distinct 
representative deformation processes. 



I. INTRODUCTION 

A. The local density of states 

Consider a system that is described by an Hamiltonian 
H(Q, P; x) where (Q, P) are canonical variables and x is 
a constant parameter. Our interest in this paper is in the 
case where (Q, P) describe the motion of a particle inside 
a cavity, and x controls the deformation of the confining 
boundary. The ID version of a cavity, also known as 
'potential well', is illustrated in Fig.l. However, we are 
mainly interested in the case of chaotic cavities in d > 1 
dimensions. Cavities in d = 2 dimensions, also known as 
billiard systems, are prototype examples in the studies of 
classical and quantum chaos, and we shall use them for 
the purpose of numerical illustrations. 
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FIG. 1. The shape of a cavity in d dimensions is defined 
by its d — 1 boundary. The confining potential is V(Q). The 
figure illustrates V(Q) for ID well. It also can be regarded 
as a cross section of d > 1 cavity. The kinetic energy of 
the particle is E — ^mv 2 . The walls of the cavity exert 
a field of force / on the bouncing particle. The hard wall 
limit corresponds to / — > oo and Vo — > oo. For theoretical 
considerations it is convenient to assume that / and Vo are 
large but finite. Mathematically it is also convenient to think 
of the embedding space as having some huge but finite volume 
(not illustrated). 



The eigenstates of the quantized Hamiltonian are 
\n(x)) and the corresponding eigen-energies are E n (x). 
The eigen-energies are assumed to be ordered, and the 
mean level spacing will be denoted by A. We are inter- 
ested in the parametric kernel 



P(n\m) = \(n(x)\m(x ))\ 2 = trace(p„p„ 



(1) 



In the equation above p m (Q,P) and p n (Q,P) are the 
Wigner functions that correspond to the eigenstates 
\m(xo)) and |n(a;)) respectively. The trace stands for 
dQdP / (2irh) d integration. The difference x — xq will be 
denoted by 5x. We assume a dense spectrum, so that 
our interest is in 'classically small' but 'quantum me- 
chanically large' energy scales. It is important to realize 
that the kernel P(n\m) has a well defined classical limit. 
The classical approximation (see remark ffl) is obtained 
by using microcanonical distributions instead of Wigner 
functions. 

Fixing n, the vector P(n\m) describes the shape of 
the n-th eigenstate in the Ho = H(Q, P; xq) representa- 
tion. By averaging over several eigenstates one obtains 
the average shape of the eigenstate (ASOE). We can also 
identify P(n\m) as the local density of states (LDOS), 
by regarding it as a function of n, where m is considered 
to be a fixed reference state. In the latter case an aver- 
age over few m-states is assumed. We shall denote the 
LDOS by P{r) where r = (n-m). The ASOE is just 
P(-r). Note that the ASOE and the LDOS are given 
by the same function. One would have to be more care- 
ful with these definitions if TL$ were integrable while Tl 
non-integrable. 

A few words are in order regarding the definition of the 
LDOS, and its importance in physical applications. The 
LDOS, also known as strength function 0-0], describes 
an energy distribution. Conventionally it is defined as 
follows: 
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LDOS(£) = (m\$tG(E)\m) 

= ^P(n\m) S(E - E n 



(2) 



where G(-E') = l/(E — H+iO) is the retarded Green func- 
tion. We are interested in chaotic systems, so it should 
be clear that our P(r) is related by trivial change of vari- 
able (E i-> r) to the above defined LDOS(.E). Our P(r) 
also incorporates an average over the reference state. The 
LDOS is important in studies of cither chaotic or com- 
plex conservative quantum systems that are encountered 
in nuclear physics as well as in atomic and molecular 
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physics. Related applications may be found in meso- 
scopic physics. Going from Ho to Ti may signify a phys- 
ical change of an external field, or switching on of a per- 
turbation, or a sudden-change of an effective-interaction 
(as in molecular dynamics 0). The so called 'line shape' 
of the LDOS is important for the understanding of the as- 
sociated dynamics. It is also important to realize that the 
LDOS is the Fourier transform of the so-called 'survival 
probability amplitude' (see j8[ for concise presentation 
of this point). 

B. Parametric Evolution 

Textbook Q formulations of perturbation theory can 
be applied in order to find the LDOS. Partial summa- 
tions of diagrams to infinite order can be used in or- 
der to get an improved Lorentzian-type approximation. 
However most textbooks do not illuminate the limita- 
tions and the subtleties which are involved in using the 
conventional perturbative schemes. It is therefore inter- 
esting to take a somewhat different approach to the study 
of LDOS. The roots of this alternate approach can be 
traced back to the work of Wigner jlo| regarding a sim- 
ple banded random matrix (BRM) model Ti = E + <5xB. 
Here E is a diagonal matrix whose elements are the or- 
dered energies {E n }, and B is a banded matrix. The 
study of this model can be motivated by the realiza- 
tion that in generic circumstances it is possible to write 
H(Q,P;x) w H (Q,P) + 5xF(Q,P). Using a simple 
semiclassical argument [|ll] it turns out that the matrix 
representation of any generic J 7 , in the eigen-basis that is 
determined by the chaotic Hamiltonian Ho, is a banded 
matrix. 

The important ingredient (from our point of view) in 
the original work by Wigner, is the emphasis on the 
parametric-evolution (PE) of the LDOS. The LDOS de- 
scribes an energy distribution: For 5x = the kernel 
P(r) is simply a Kroneker delta function. As 5x be- 
comes larger, the width as well as the whole profile of 
this distribution 'evolves'. Wigner has realized that for 
his WBR model there are three parametric regimes. For 
very small Sx we have the standard perturbative struc- 
ture where most of the probability is concentrated in 
r = 0. For larger Sx we have a Lorentzian line shape. 
But this Lorentzian line shape does not persist if we fur- 
ther increase Sx. Instead we get a semi-circle line shape. 
Many works about the LDOS have followed [0-[|, but the 
issue of PE has not been further discussed there. The em- 
phasis in those works is mainly on the case where Ho is 
an intcgrablc or non-interacting system, while H is pos- 
sibly (but not necessarily) chaotic due to some added 
perturbation term. 

The line of study which is pursued in the present work 
has been originated and motivated by studies of quan- 
tum dissipation [p^|-pT[. Understanding PE can be re- 
garded as a preliminary stage in the analysis of the energy 



spreading process in driven mesoscopic systems. Note 
that the LDOS gives the energy re-distribution due to a 
'sudden' (very fast) change of the Hamiltonian. Unlike 
the common approach for studies of LDOS, we assume 
both H and Ho to be chaotic. Both correspond to the 
same parametrically dependent Hamiltonian H(Q, P; x), 
and there is nothing special in choosing a particular value 
x = xq as a starting point for the PE analysis. 

C. Main results 

The theory of PE, as discussed in general in |L^-|l4|] 
and in particular in ||,[l5| takes us beyond the random- 
matrix-theory considerations of Wigner. There appear 
five (rather than two) different parametric scales (see re- 
mark pj). These are summarized by Table 1. 



Sxf Is it possible to linearize TL(x + Sx)? 

SXc" 1 Is it possible to use standard perturbation theory? 

Sx vrt Do perturbative tail regions survive? 

5a; N u Do non-universal core features show up? 

5x S c Is it possible to use semiclassical approximation? 



TABLE I. The parametric scales in the general the- 
ory of PE are listed (left column) along with the ques- 
tions that motivate their introduction. The distribution P(r) 
may contain perturbative tail regions (for Sx <C Sx pTt ), and 
non-perturbative core regions (for Sx > Sx^" 1 ). Non-universal 
(system specific) features may manifest themselves in the core 
structure for Sx S> Sx NU . In generic examples Sx ^> Sx sc 
allows classical approximation for P(r). We are going to ex- 
plain that only two independent parametric scales survive in 
the hard wall limit. 

In the present paper we consider cavities with hard 
walls. We are going to explain that because of the 'hard 
wall limit' there are only two independent parametric 
scales: One is Sx^" 1 and the others (see remark ||) coin- 
cide with Sx NV . Assuming that Sx^™ 1 and S well 
separated, it follows that there are three distinct para- 
metric regimes in the PE of our system. These are the 
standard perturbative regime (Sx <C Sxl'"), the core- 
tail regime (Sx^." 1 <C Sx <C Sx m ), and the non- universal 
regime (Sx Sx m ). 

The exploration of the three parametric regimes in the 
PE of a deformed cavity with hard walls is the main issue 
of the present paper. To the best of our knowledge such 
detailed exploration has not been practical in the past. 
We owe our ability to carry out this task to a new pow- 
erful technique for finding clusters of billiard eigenstates 
|]l9| , f20f . There are also some secondary issues that we are 
going to address: 

(a) In the strict limit of hard walls the PE becomes 
non-unitarity. We shall use the ID well example in order 
to shed light on this confusing issue. In particular we 
demonstrate that any truncation of the PE equation leads 
to false unitarity due to a finite-size edge effect. 
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(b) For special deformations, namely those that consti- 
tute linear combination of translations rotations and dila- 
tions, the parametric scales Sx*™ 1 and S x N u coincide. Con- 
sequently there is no longer distinct core-tail regime, and 
the PE features a quite sharp transition from the stan- 
dard perturbative regime to the non-universal regime. 

(c) In the non-universal regime we demonstrate that 
our numerical results are in accordance with our theo- 
retical expectation ||. Namely, the width of the LDOS 
profile is determined by time-domain semiclassical con- 
siderations, rather then by phase-space or random-wave 
considerations. 

(d) The last section puts our specific study in a larger 
context. We explain why Wigner's scenario of PE is not 
followed once hard walls are considered. 



II. THE CAVITY SYSTEM 

We consider a particle moving inside d dimensional 
cavity whose volume is V. The kinetic energy of the 
particle is E — |mt> 2 , where m is its mass, and v is 
its velocity. It is assumed that this motion is classically 
chaotic. The ballistic mean free path is l hl . One can use 
the estimate £ bl ~ V/A, where A is the total area of the 
walls. The associated time scale is r bl = £ bl /v. 

The penetration distance upon a collision is £ = E/ f, 
where / is the force that is exerted by the wall. Upon 
quantization we have an additional length scale, which 
is the De-Broglie wavelength A B = 2-K%/(mv). We shall 
distinguish between the hard walls case where we assume 
I < A B *C £bu an d soft walls for which A B <C £■ Note that 
taking h — > implies soft walls. 

There is a class of special deformations that are shape- 
preserving. These are generated by translations, rota- 
tions and dilations of the cavity. A general deformation 
need not preserve the billiard shape nor its volume. We 
can specify any deformation by a function -D(s), where s 
specifies the location of a wall element on the boundary 
(surface) of the cavity, and D(s)8x is the normal dis- 
placement of this wall element. In many practical cases 
it is possible to use the convention |D(s)| ~ 1. With this 
convention Sx has units of length, and its value has the 
meaning of typical wall displacement. 

The eigen-energies of a particle inside the cavity are 
in general x-dcpcndcnt, and can be written as E n — 
(Hk n ) 2 / (2m) . The mean level spacing is 



(3) 



where Qd — 2n, Ait, ■ ■ ■ for d = 2, 3, • • •. In our numerical 
study we shall consider a quarter stadium with curved 
edge of radius 1 and straight edge of length 1. The 
'volume' of the quarter stadium is V = 1 + 7r/4. The 
'area' of its boundary A = 4 + ir/2 is just the perimeter. 
We shall look on the parametric evolution of cigenstates 
around k ~ 400 where the mean level spacing in k units 
is A = A/(h,v) re 0.0088. 



III. PARAMETRIC EVOLUTION 

Consider the quantum-mechanical state tj) = \m(xo)}. 
We can write ip — En a n(x)\n(x)) . The parametric ker- 
nel can be written as P(n\m) = \a n (x)\ 2 . It is a standard 
exercise to obtain (from dip j dx = and differentiating by 
parts) the following equation for the amplitudes: 



da n 
dx 



(4) 



In order to get P(n\m) one should solve this equation 
with the initial conditions a n (xo) — 5 nm . The transitions 
between levels are induced by the matrix elements 



W 



ih 



E n —E„ 



dx 



(5) 



and we use the 'gauge' convention W„ m =0 for n=m. 
(Only one parameter is being changed and therefore 
Berry's phase is not an issue). 

Eq.([|) is a possible starting point for constructing a 
perturbation theory for the PE of P(r). See Ref. (IJ] for 
more details. As an input for this equation we need the 
matrix elements of dTL/dx. These can be calculated using 
a simple boundary integral formula Jl6[ whose simplest 
derivation [Q is as follows: The position of the particle 
in the vicinity of a wall element is Q — (z,s), where s 
is a surface coordinate and z is a perpendicular 'radial' 
coordinate. We take / = oo so that 



&H 

dx 



= -D(s) V 6(z) 



(6) 



The logarithmic derivative of the wavefunction on the 
boundary is ip(s)/tp(s) where (p(s) = n-Vip, and n is a 
unit vector in the z direction. For z > the wavefunc- 
tion tp(Q) is a decaying exponential. If Vq is large enough, 
then the exponential decay is fast, and we can treat the 
boundary as if it were locally flat. It follows that the 
logarithmic derivative of the wavefunction on the bound- 
ary should be equal to — v / 2mVb/?i. Consequently one 
obtains the following expression for the matrix elements: 



m 

dx 



2m 



j> ip n (s)tp m (s) D(s)ds 



(7) 



In the one-dimensional case the boundary integral is re- 
place by the sum J2 S L Pn(s)(p m (s)D(s) where s = 1, 2 are 
the two turning points of the potential well. 



IV. HARD WALLS AND NON-UNITARITY 

For the purpose of the following argumentation it is 
convenient to take / = oo, but to keep Vq large but 
finite. Mathematically it is also convenient to think of the 
embedding space as having some huge but finite volume. 
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We would like to illuminate a subtlety which is associated 
with the hard wall limit Vq — > oo. For any finite Vo the 
parametric kernel satisfies 



^P(n\m) = p t , 



(8) 



is supported by the simple pedagogical example of the 
next section. 



V. PARAMETRIC EVOLUTION FOR A ID BOX 



withp total = 1. This follows from the fact that \n(x)) is a 
complete orthonormal basis for any x. However, for hard 
walls (Vq = oo) this statement is not true. This implies 
that for hard walls the PE in non-unitary. We are going 
to explain this point below. 

Let us denote the volume of the original cavity by Vo 
and of the deformed cavity by V. The volume shared by 
the deformed and the undeformed cavities will be denoted 
by V n V and we shall use the notation i] = (V n V) / V . 
For the purpose of the following argumentation let us 
consider a reference state m whose energy E m is well be- 
low Vq. Let us also assume that the wall displacement is 
large compared to De-Broglie wavelength. Consequently 
the expression for P(n\m) has the following semiclassical 
structure: 



Consider a ID box with hard walls, where the free mo- 
tion of the particle is within < Q < a. The eigenstates 
of the Hamiltonian are 



*(a)> 



{-l) n \ - sin(fc„Q) 



(10) 



where k n = nx (n/a) is the wavenumber, and n — 1,2, ... 
is the level index. The phase factor (—1)" has been intro- 
duced for convenience. We consider now the parametric 
evolution as a function of a. One easily obtains 



(n(a)|m(a )) = (-l) n 0? 



sin(7r?7n) 



2m 



2 2 



(11) 



where i] = a /a is assumed to be smaller than 1, corre- 
sponding to expansion of the box. The probability kernel 
P(n\m) w X] x f(E n -E m ) + (I-77) x g(E n -E c ) (9) is p( n | m ) = |( n |m)| 2 . One can verify that the parametric 



where E c — (Vo + E m ) ~ Vo. The above result can 
be deduced by assuming that the wavefunctions look er- 
godic in space, but still that they are characterized by a 
well-defined local wavelength. An equivalent derivation 
is obtained by using the phase-space picture of Q . Both 
/ and g in the above expression have unit normalization, 
and therefore p total — 1 for any finite Vo- However, for 
hard walls (Vq = 00) we have E c — 00 and therefore 
Ptotai = V- We may say the the operation of taking the 
hard wall limit does not commute with the summation 
in Eq.(||). An analogous statement can be derived re- 
garding the summation ^2 m P(n\m), with the respective 
definition 7/ = (V n V) /V. 

The correctness of the above observation becomes less 
trivial if we consider (Q) with expressions (||) and (Q) sub- 
stituted for the matrix elements. Looking on (||) with ([?]) 
it looks as if the matrix W nm is hermitian, and therefore 
should generate unitary PE. But this statement is math- 
ematically correct only for (any) finite truncation N of 
the PE equation. For N = 00 the matrix W„ m becomes 
non-hermitian. It turns out that for any finite N, there 
is a pile-up of probability in the edges of the spreading 
profile, due to finite-size effect. We shall demonstrate 
this effect in the next section using a simple ID example. 
In other words, if we solve Eq.(Q) for hard- walled cavity, 
we get as a result Eq.(|]) with E c = En- For N — 00 we 
get E c = 00 and therefore p to tai = 77 in accordance with 
the conclusion of the previous paragraph. 

Thus if either Vo < 00 or N < 00 then we have unitary 
PE. But for hard walls, meaning Vo — 00 with N = 00, we 
have non-unitary PE. The lost probability is associated 
with the second term in Eq.(Bh. This term is peaked 
around a high energy E c . For hard walls E c — 00 and 
consequently some probability is lost. The above picture 



evolution in the ao i— > a direction is unitary, meaning that 
J2n P{n\m) = 1. On the other hand in the a n ao di- 
rection the parametric evolution is non-unitary, because 
P( n \ m ) — V- The profile of P(n\m) for fixed n is 
illustrated by a dashed line in Fig. 2. 

We can restore unitarity by making Vq large but finite. 
In such variation of the above calculation leads 

to the following picture: Consider the overlap of a ref- 
erence level n(a) with the levels rn(ao). As in the case 
Vo = 00 there is a probability 77 which is located in the 
levels whose energies are E m w E n . But now the "lost" 
probability (1 — 77) is located in the levels whose energies 
are E m w E n + Vo. Thus we have J2 m P(n\m) = 1 rather 
than J2 m P( n \ m ) = V- 

We consider again the case Vo = 00. The normal 
derivative on the boundary is if n (a) = \J (2/a) k n . Hence 
we can easily get the following result 



I Wnm = 



1 2nm 

■^n(a)<p m (a) = -7 (12) 



It is more convenient to use a = ln(a) for parameteriza- 
tion. Hence the equation that describes the parametric 
evolution is 



da n ^ 



2nm 



(13) 



For any finite truncation N < 00 this equation manifestly 
generates unitary parametric evolution. It is only for 
N = 00 that it becomes equivalent to the non-unitary 
evolution of the ID box. Again, one can wonder where 
the 'lost' probability is located if A < 00. The answer 
is illustrated in Fig. 2. We see that the 'lost' probability 
piles up at the edge of the (truncated) tail. 
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VI. MATRIX ELEMENTS FOR CHAOTIC 
CAVITY 




FIG. 2. Upper subfigure: An image of the kernel P(n\m) 
for 13% expansion of the box (ie a/ao = 1.13). The ker- 
nel has been calculated numerically using Eq.(h3J) with finite 
truncation N = 256. Lower subfigure: The profile of a rep- 
resentative row of P(n\m). The dashed line is the N — oo 
analytical result using Eq.(|ll|). 



It is possible to use semiclassical considerations O] 
in order to determine the band profile of the matrix 
Eq. ((?]). The application to the cavity example has been 
introduced in and numerically demonstrated in [ jl7| . 
The accuracy of this semiclassical estimate is remarkable. 
Here we summarize the recipe. First one should gener- 
ate a very long (ergodic) trajectory, and define for it the 
fluctuating quantity 



OH 

dx 



^2m« cos(6» col ) D col 5{t-t c , 



(14) 



where t col is the time of a collision, D coi stands for D(s) 
at the point of the collision, and i>cos(# col ) is the normal 
component of the particle's velocity. Each delta spike (for 
soft walls it is actually a narrow rectangular spike) cor- 
responds to one collision. Now one can calculate the cor- 
relation function C(t) of the fluctuating quantity Tit), 
and its Fourier transform C(lo). The semiclassical esti- 
mate for the band profile is 



&H 

dx 



A ^ I E n —E„ 



(15) 



Ref. jl8[ contains a systematic study of the function 
C{u>). For large lo, meaning lo 3> 1/T b] , one can use 



2mV(|cos0| 3 )y (j)[D(s)} 2 ds I K.) 



where the geometric factor is (| cos(6 , )| 3 ) = l,4/37r, ••• 
for d = 1, 2, • • •. A lengthy calculation |14[ reveals that 
Eq. (p~5|) with ( |l6| ) substituted, is an exact global result if 
we could assume that the cavity eigenfunctions look like 
'random waves', and that different wavefunctions are un- 
corrected. However, it turns out that to take this 'ran- 
dom wave' result as a global approximation is an over- 
simplification. For lo <C 1/tm, using the semiclassical 
recipe and assuming strongly chaotic cavity, one obtains 



C(w) « C*(oo) x (r bl ujy 



(17) 



with 7 = 4 for dilations and translations, 7 = 2 for ro- 
tations, and 7 = for normal deformations. We use the 
term 'special deformations' 18 in order to distinguish 
those deformations that has the property C(u>) — > in 
the limit lo — > 0. Any combination of dilations, trans- 
lations and rotations is a special deformation. Around 
the bouncing frequency (lo ~ l/ 7 ") the function C(lo) 
typically displays some non-universal (system and defor- 
mation specific) structure. This is true for any typical de- 
formation, but for some deformations the non-universal 
features are more pronounced. If the cavity has bouncing 
ball modes, we may get also a modified (non-universal) 
behavior in the small frequency limit. For the purpose 
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FIG. 3. Band profiles for deformations of the quarter sta- 
dium as denned by Eq.(^|). Upper plot: Rotation around 
the stadium center. Lower plot: Generic (non-special) de- 
formation involving displacement of the curved edge. It is 
important to notice that for the special deformation we have 
C{oj) — > in the limit ui — * 0. 



of general discussion it is convenient to assume that the 
interpolation between (17) and ( |l6| ) is smooth, but in ac- 
tual numerical calculation the actual C(lv) is computed 
(see below). 

As a numerical example we have picked the sta- 
dium billiard. We have found the eigenstates of a de- 
symmetrized (quarter) stadium as described in Ref. [p"7| . 
We have selected those eigenstates whose eigen-energies 
k n are in the vicinity of k — 400. Our two representative 
deformations are: (a) rotation around the stadium cen- 
ter; and (b) generic (non-special) deformation involving 
the curved edge. In the latter case the curved edge of the 
quarter stadium (0 < s < tt/2) is pushed outwards with 
D(s) — (cos(s)) 2 , while for the straight edges D(s) = 0. 
(The corner s — is the 90° intersection of the curved 



edge with the long straight edge). The respective band 
profiles are displayed in Fig. 3. The band profile has been 
defined as 



ip„(s)<p m (s) D(s)ds 



(18) 



where k — (k n — k m ) is the distance from the diagonal. 
Note that B(n) is just a scaled version of the semiclassical 
C(ui) as implied by Eq.(|l5|) with (fjj). The remarkable 
agreement of B(n) with the semiclassical calculation has 
been demonstrated in [^JTsJ. 

It is important to realize that in the hard wall limit 
(which is assumed here) the matrix (dTLj dx) nm is not 
a banded matrix. It would become banded if we were 
assuming soft walls. For soft walls C(ui) becomes vanish- 
ingly small for u> 3> v/£. The bandwidth in energy units 
is Af, = hv/i, and in dimcnsionlcss units it is 



b = =^ 



A 



V 



I A 



d-l 



(19) 



Unless stated otherwise we have b = oo. 



VII. PARAMETRIC EVOLUTION 
RESULTS 



NUMERICAL 



The parametric evolution of P{r) for rotation and for 
generic deformation of the stadium is illustrated by the 
images of Fig. 4 and by the plots of Figs. 5-6. The calcula- 
tion of each P(r) profile is carried out as follows: Given 
Sx we use the method which is described in pcfl in order 
to calculate the matrix P(n\m). Then we plot the ele- 
ments of P(n\m) versus k — (k n (x) — k m (0)). In order 
to obtain the average profile theplot is smeared using 
standard procedure (see remark |2l]]). The transforma- 
tion from k to T = [n — m) is done using the relation (see 
remark f2lf): 



A- 



1, 5V 

d kx V 



(20) 



Above A is the mean level spacing of the {k n } spectrum, 
and <5V is the volume change that is associated with the 
deformation (it is approximately proportional to Sx). If 
the deformation is volume preserving (as in the case of 
rotation) then the second term equals zero. But for the 
generic deformation that we have picked in our second 
numerical example, the volume is not preserved, and the 
systematic 'downwards' shift of the levels should be taken 
into account. 

Looking first in the case of rotation, we see clearly two 
parametric regimes: The standard perturbative regime 
(Sx < 0.2), and the non-universal regime (6x > 0.2). Let 
us clarify this observation. We see that for Sx < 0.2 most 
of the probability is well concentrated in r = 0. This 
implies that we can use standard perturbation theory in 
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FIG. 4. Each column is an image of P(r) versus k for a 
different value of 8x. There are 41 columns. The value Sx ~ 1 
corresponds roughly to A B displacement of the boundary. The 
upper figure is for rotation, and the lower is for generic de- 
formation. Note that the r — component is excluded from 
the image. Instead we have plotted over the image an r = 
line wherever this component contains more than 50% of the 
probability. In the lower figure this line cannot be resolved 
from the developing core region. 



order to estimate the r ^ probabilities. On the other 
hand for Sx > 0.2 the perturbative nature of P(r) is 
destroyed. Now P(r) becomes smoother, and eventually 
(for Sx > 0.5) there is a very good fitting with Lorentzian 
(see lower plot in Fig.5). 

The qualitative explanation for the Lorentzian profile 
is as follows: For Sx > 0.5 the typical displacement of 
the walls is of the order of A B . Therefore the \n(x)) 
eigenstates become uncorrelated with the |m(0)) eigen- 
states. Consequently P(r) becomes Sx independent. The 
Lorentzian profile agrees with the assumption of uncor- 
related random waves as explained in Appendix A. 

Let us look now in the case of generic deformation. 
Here we see clearly three parametric regimes: The stan- 
dard perturbative regime (Sx < 0.004), the core-tail 
regime (0.004 < Sx < 0.2), and the non-universal regime 
(Sx > 0.2). Let us clarify this observation. As in the 
case of rotation there is a standard perturbative regime 
(Sx < 0.004) where most of the probability is well con- 
centrated in r — 0. For larger deformation, namely for 
Sx > 0.004, standard perturbation theory is no longer 
applicable because the n = m level is mixed non per- 
turbatively with other (neighboring) levels. As a re- 
sult P(r) contains a non-perturbative 'core' component. 
However, for 0.004 < Sx < 0.2 we definitely do not get a 
Lorentzian. Rather the tails of P(r) keep growing in the 
same way as in the standard perturbative regime. 

In the case of the generic deformation, as in the case 
of rotation, we enter the non-universal regime, and even- 
tually (for Sx > 0.6) we get a smooth Lorentzian-like 
distribution. However, the Lorentzian-like distribution is 
not identical with that of the rotation case. Also the sim- 
ilarity to proper Lorentzian is far from being satisfactory, 
(see lower plot in Fig. 6). This means that the random- 
wave picture of Appendix A is an oversimplification. 

In the following sections we are going to summarize the 
theoretical considerations || that explain the observed 
parametric scenario. In particular we are going to illumi- 
nate the way in which non-perturbative features emerge; 
to clarify the crossover to the non-universal regime; and 
to explain the specific nature of the non-universal distri- 
bution. 



VIII. THE STANDARD PERTURBATIVE 
REGIME 

Standard perturbation theory gives the following first 
order expression for the LDOS 

P(n\m) « S nm + Sx 2 \(dn/dx) nm \ 2 

(£j n — £j m ) 

This expression is most straightforwardly obtained by in- 
specting Eqs.Q and (|B|). We can define the (total) tran- 
sition probability as 

p(Sx) = £p(r) (22) 
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FIG. 5. Representative plots of P(r) for the case of rota- 
tion. The upper subfigure is for 0.0010 < Sx < 0.2512. The 
r = component is excluded. In the lower subfigure the fitted 
Lorentzian is indistinguishable from the actual profile. 
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FIG. 6. Representative plots of P(r) for the case 
of generic deformation. The upper subfigure is for 
0.0010 < Sx < 0.2512. In the lower subfigure a fitted 
Lorentzian is overlayed for the purpose of comparison. 
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Using Eq.(|2 
estimate: 

p{5x) 



combined with (|15[) we get the following 



Sx 2 



du> C (w) 
2^ up- 



(23) 



Standard perturbation theory is applicable as long as 
p(Sx) <C 1. This can be converted into an equivalent 
inequality Sx <C Sx^" 1 . By this definition Sx^" 1 is the 
parametric deformation which is needed in order to mix 
the initial level m with other levels n m. 

If we use Eq.(23) for a special deformation, then we 
have 7 > 1, and consequently the integral is not sensitive 
to the exclusion of the \r\ < A region. As a result we 



have Sx*™ 1 cx h. Using Eq.Q with we get 



8x« 



I special 



[D{s)} 2 ds 



x A F 



(24) 



In case of generic deformation Eq.d24) is not valid be- 
cause the value of the integral in (p3[) is predominantly 
determined by the to ~ A/Ti lower cutoff rather than by 
an effective u> ~ l/r bl lower cutoff. As a result one ob- 
tains 



Sx* 



generic 



£ H i 0[D(s)]-,/s 



x A F 



(25) 



where £ H = vt H is the length which is associated with the 
Heisenberg time t H — 2nT%/ A. 

It is illuminating to use the convention D(s) ~ 1, such 
that Sx measures the typical displacement of a wall el- 
ement. With this convention we get from Eqs.(|2~4|)-(l2T 
the following: 



Sxl 



Sx 



I special 



A P 



generic 
^ x c Idiffractive 



\d-l 



A efl 



x A F 



Ap 



(26) 

(27) 
(28) 



In the generic case the effective area of the deformed 
boundary A eff may be smaller than the total area A of 
the boundary. The effective area A cff can be formally de- 
fined by comparing Eq.@ with Eq.©. Eq.(||) has 
been added for sake of completeness of our presentation. 
It correspond to the diffractivc limit A off — > X^ 1 . Note 
that this limit is beyond the scope of the present study. 
Thus in the generic case, the wall displacement needed 
to mix levels is much smaller than A B - In the generic 
case Sxl™ 1 cx ^( d+1 )/ 2 rather than Sx^ m cx h. What hap- 
pens with perturbation theory beyond 5x^1 This is the 
subject of the next section. 



IX. THE CORE- TAIL REGIME 

For 5x > Sx^. 1 " standard perturbation theory diverges 
due to the non-perturbative mixing of neighboring levels 



on small scale: Once 5x becomes of the order of Sx^" 1 
several levels are mixed, and as Sx becomes larger, more 
levels are being mixed non-perturbatively. Consequently 
it is natural to distinguish between core and tail regions 
]T3| fl5[ . Most of the spreading probability is contained 
within the core region, which implies a natural extension 
of first-order perturbation theory (FOPT): The first step 
is to transform Eq. (Q) into a new basis where transitions 
within the core are eliminated; The second step is to use 
FOPT (in the new basis) in order to analyze the core- 
to-tail transitions. Details of this procedure, which is 
in the spirit of degenerate perturbation theory, were dis- 
cussed in Jl4| . The most important (and non-trivial) con- 
sequence of this procedure is the observation that mixing 
on small scales does not affect the transitions on large- 
scales. Therefore we have in the tail region P(n\m) cx Sx 2 
rather than P(n\m) cx Sx. The validity of this observa- 
tion has been numerically illustrated in Jl5| . 

Following the above reasoning we define the tail region 
as consists of those levels whose 'occupation' can be cal- 
culated using perturbation theory, while the core is the 
non-perturbative component in the vicinity of r = 0. As- 
suming that only one scale characterize the core width, 
one arrives at the following practical approximation: 



En— En 



Sx 2 



{T{5x)Y + (E n -E m y 



(29) 



It is implicit in this definition that (E n —E m ) should be 
regarded as a function of r = (n — m). The parameter 
T(5x) is determined (for a given Sx) such that P prt (r) has 
a unit normalization. One may say that T(5x) regularizes 
the behavior around r = 0. For generic deformation we 
get 



T(Sx) = b A 



( Sx 



\Sxl 



x A 



(30) 



The core is defined as the region \r\ < bo, and the outer 
(|r| > bo) regions are the tails. For Sx <C <Sx^ m we 
get 6 « 1 and the core-tail structure Eq.(p^) becomes 
equivalent to the standard perturbative result Eq. ( ^l| ) . It 
should be clear that the core-tail structure is a general- 
ization of Wigner's Lorentzian. It is indeed a Lorentzian 
in the special case of a 'flat' band profile. 

We turn now to analyze our numerical results. We 
have verified (see eg Fig. 7) that for Sx < 0.2 we have 
good agreement with perturbation theory irrespective of 
whether we have a core component (which is the case for 
the generic deformation) or not (which is the case for the 
rotation). As we come closer to Sx ~ 0.2 the agreement 
becomes worse, and for Sx > 0.3 we have a total col- 
lapse of perturbation theory. In Fig. 8 we display the to- 
tal transition probability p(Sx) as a function of Sx. This 
plot should be used in order to numerically determine the 
value of Sx^,™, say as the value where p(5x) = 1/2. 

Fig. 8 also displays comparison with the cor resp onding 
perturbative calculation (using Eq.(|22|) with(|29|)). The 
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FIG. 7. Comparison between the perturbative calculation 
fprt(r) and the actual P(r). The upper subfigure is for the 
rotation and the lower subfigures are for the generic deforma- 
tion. The lower subfigure gives a zoom over the the 'birth' of 
the non-perturbative semiclassical saturation profile. 




FIG. 8. The upper and the middle subfigures are p(Sx) 
for the rotation and for the generic deformation respectively. 
The lower subfigure is the perturbative calculation of the core 
width bo for the generic deformation. Also displayed is the 
calculated b c for P prt (r), and the actual b c for P(r). 
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agreement (for Sx < Sx^°) in case of the rotation is re- 
markable. For larger values of Sx one observes a linear 
drop of the total probability due to the non-unitary na- 
ture of the evolution. (This drop is clearly linear in a 
linear- linear scale which is not displayed). The eventual 
rise of the total transition probability (beyond 1) in case 
of the generic deformation, reflects numerical errors in 
the determination of the small k overlaps. See J2(J for 
further details. 

The lower subfigure in Fig. 8 displays the calculated 
core width b for the generic deformation. Recall that 
b is determined, given Sx, such that P PTt (r) of Eq.(gg) 
is normalized. Also calculated is the width b c of the re- 
gion that contains 50% of the probability. The width 
b c is calculated for both P VTt (r) and P(r). Note that 
b c = 1 in the standard perturbative regime. The de- 
termination of Sx"}" 1 for the generic deformation becomes 
more convenient by using this plot. Also the crossover 
(at 5x NU ~ 0.2) to the non-universal regime is most pro- 
nounced. 



X. THE NON-UNIVERSAL REGIME 

The validity of Eq. (^9|) as a global approximation relies 
on the assumption that the core is characterized only by 
one scale (which is bo). But this assumption ceases to 
be true if Sx is large enough. We have explained in Ref. 
H that the width of the core defines a 'window' through 
which we can view the 'landscape' of the semiclassical 
analysis. As Sx becomes larger, this 'window' becomes 
wider, and eventually some of semiclassical structure is 
exposed. This is marked by the non-universal parametric 
scale Sx NV . For Sx larger than Sx NU , the non- universal 
structure of the core is exposed. 

What is the semiclassical structure of the LDOS? 
Time-domain semiclassical considerations (See Eq.(10) 
of H and related discussion there) imply that the non- 
universal structure of the core is 



P{n — m) 



1 



7T 6E% C 



6E* 



(E n —E n 



(31) 



with SE SC — h/r col . The definition of the collision rate 
1/t co1 is similar to that of l/r bI . The former is the col- 
lision rate with the deformed area A dfr , and therefore it 
may be smaller than l/r bl , because t u is related to the 
total area A. 

If, by mistake, we identified t co1 with r bl , then we would 
get the random wave result which is derived in Appendix 
A. This would be an over simplification. If it were true, 
it would imply that for Sx > Sx m we should have got the 
same Lorentzian-distribution in both cases (the rotation 
and the generic deformation). What we see, as a matter 
of fact, is that for the rotation (Fig. 5) we have a rea- 
sonably good agreement with Lorentzian whose width is 
Sk sc — 1.5, whereas for the generic deformation (Fig. 6) 
there is rough agreement with Lorentzian whose width is 



Sk sc = 0.26. The smaller width in the latter case clearly 
reflects having larger r col . For the generic deformation we 
also have pronounced non-Lorentzian features. Actually 
the global fitting to Lorentzian is quite bad. Our under- 
standing is that these features are due to the bouncing- 
ball trapping: It leads to non-exponential decay of the 
time-dependent survival probability (see Q for defini- 
tion of the latter term), and hence to the observed non 
Lorentzian features of the spreading profile. 

We turn now to explain how fa NU is determined. By 
definition it is the deformation which is required in order 
to expose features of the semiclassical landscape. These 
features start to be exposed once T(5x) ~ SE SC which 
should be converted into an equivalent expression Sx ~ 
5x NU . Thus we get 



5x K 



Adfr 

A c ff 



x Ap 



(32) 



Here A dfr is the geometric area of the deformation (in the 
sense of scattering cross-section), while A off is the effective 
area of the deformation. The definition of the latter is 
implied by comparing Eq. ( p7| ) with Eq. (^5|) . By rescaling 
-D(s) i— ► aD(s) and Sx i— ► Sx/a, we can make A olf = A dfr 
by convention. For special deformations Sx NV coincides 
with Sx^" implying that we get into the non-universal 
regime as soon as we have a breakdown of standard per- 
turbation theory. 

Our theoretical consideration so far do not imply a 
total collapse of perturbation theory. We may have in 
principle a co-existence of non-universal core component 
and perturbative tails. Actually we see such co-existence 
in the lower subfigures of Fig. 7, mainly for Sx = 0.2512. 
The right peak around k = is clearly non-perturbative, 
while the rest of the profile is in reasonable (though not 
very good) agreement with the perturbative calculation. 
We are going to explain in the next section that the total 
collapse of perturbation theory for Sx > 0.3 is actually 
not related at all to the appearance of non-universal fea- 
tures in the core structure. It is only circumstantial that 
in the hard wall limit this collapse happens as soon as we 
enter the non-universal regime. 



XI. THE COLLAPSE OF PERTURBATION 
THEORY 

A good starting for the following discussion is to con- 
sider the classical approximation for P(r). Namely, 
in Eq.(|^) one approximates p„(Q,P) and p m (Q,P) by 
microcanonical distributions. A phase-space illustra- 
tion of the energy surfaces which support p n (Q,P) and 
Pm(Q, P) can be found in Fig.l of Ref. The classical 
P(n\m) equals to the overlap of these surfaces. 

If we were dealing with a generic system we 
could introduce a linearized version of the Hamilto- 
nian H(Q, P; x) = H (Q, P) + 5xF(Q, P; Sx=0), where 
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T(Q,P\x) = dTi/dx. By definition this linearization is 
a good approximation provided 5x <C Sx^, where Sx^ is 
the classical correlation scale of .F(Q, P; x) with respect 
to x. In the classical linear regime the classical P(r) has 
the scaling property P(r) — l/5xP(r/5x). 

An equivalent definition of Sx",! in the quantum- 
mechanical case is obtained by looking on the x de- 
pendence of the matrix elements of J-(Q, P;x) in some 
fixed basis. Again we define Sx^ as the respective cor- 
relation scale. It is quite clear that for cavity with soft 
walls we have 



Sx",!; = i [soft walls] 



(33) 



where I = E/f has been defined as the penetration dis- 
tance upon collision. From purely classical point of view 
the hard wall limit i — > is a non-linear limit. But this 
is not true quantum-mechanically. Here we have 



5x1 



A B [hard walls] 



(34) 



for I < A B . The terminology 'classical correlation scale' 
while referring to Sx^ becomes misleading here, but we 
shall keep using it anyway. 

The theory of the core-tail structure jl4j is valid only in 
the linear regime Sx <C Sx^. Let us assume for a moment 
soft walls. We can ask: Is Sx <C fo^ 1 a sufficient condition 
for having a core-tail structure? The answer is definitely 
not. Perturbation theory has a final collapse once the 
core width bo(8x) becomes of the order of the band width 
b. This defines a parametric scale Sx plt . For 8x > 5x plt 
the LDOS becomes purely non-perturbative. In Wigner's 
theory of random banded matrices this corresponds to 
the crossover from Lorentzian to semicircle line shape. 

However, in the limit of hard walls the above mech- 
anism of collapse becomes irrelevant because the band 
width is infinite (b = oo). On the other hand we still have 
to satisfy the inequality Sx <C Sxf. Thus, for Sx > A B we 
expect a total collapse of perturbation theory, as indeed 
observed in the numerical study. 



APPENDIX A: OVERLAP OF UNCORRELATED 
RANDOM WAVES 

It is possible to estimate the overlap |(n|m)| 2 if 
we assume that |n) and \m) are uncorrelated random- 
superpositions of plane-waves: A random-superpositions 
of plane-waves is characterized by the correlation func- 
tion 

(V>r(xi)^r(x 2 )) = ~ Cos(fc|x 2 - Xi|) 

where Cos(fcr) = (exp(ik • r))n is a generalized Bcsscl 
function (for further details see Appendix D of jbi]]). As- 
suming that the wavefunction of |n) is uncorrelated with 
the wavefunction of |m) one obtains 



<IH™}| 2 ) = [y 



J J Cos(fc„|x 2 — xi|) Cos(fc TO |x 2 — xi|) dx 1 dx2 



The integration is over the whole volume of the cavity. 
Using the definition of the Cos function we can cast this 
expression into the form (|(ro|m)| 2 ) = (/(q)) q where the 
average is over the difference q = (k 2 — ki), with all 
possible orientations for |ki| = fc m and for |k 2 | = fc n . 
The function /(q) is defined as follows: 



/(q) 



e iqx dx 



The function /(q) depends mainly on q = |q|. We can 
obtain an estimate for f(q) by considering a spherical 
cavity |x| < R with the same volume. Using spherical 
coordinates it is straightforward to obtain the following: 



/(<*) 



TT 



1 



d-l 



(R 2 



2N(d-l)/ 2 



above J-T is a Fourier transform from x to q. For q <C 
1/R we have simply f(q) — (V) 2 . For q ^> 1/R we have 
f(q) = (V) 2 / (Rq) d+1 . This is because of the singularity 
at x — ±R. The average over q can be done using again 
spherical coordinate. We have q « ((Afc) 2 + 2fc 2 (l — 
cos#)) 1/2 , where Afc = |fc n — k m \, and we can transform 
the dd integration into dq integration: 



</(q)) c 



d-2 



d-l 




(Afc) 2 )((2fc) 2 



d-3 



2fc 2 



where is the solid angle in d dimensions. There is no 
point in trying to carry an exact integration. Rather, it 
is important to observe that a practical approximation 
for the overlap is: 



(|(n|m)| 5 



1 

kR 



d-l 



1 



l + (i?-(fc„-fc m )) 2 



In the latter expression we have neglected the d- 
dependent normalization prefactor. It is a 'practical' ap- 
proximation since it gives the correct behavior for both 
small and large values of Afc. The interpolation around 
Afc ~ 1/R cannot be trusted. 
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The classical approximation for P(n\m) works well for 
ergodic eigenstates. As explained later in the text we ac- 
tually consider the averaged profile. Namely, P(n\m) is 
regarded as a function of r = n — m, and it is aver- 
aged over the reference state m. With this procedure the 
effect of the minority of non-ergodic eigenstates can be 
neglected in the classical limit. 

The two parametric scales of Wigner's RMT model cor- 
respond to SXc™ and 8x plt of Table 1. Consequently there 
are three parametric regimes in Wigner's theory. These 
are the standard perturbative regime (Sx <C Sx^" 1 ), the 
Lorentzian regime (Sx^" 1 <C 8x <C Sx plt ), and the semi- 
circle regime (<5:r S> 8x pT t)- As an artifact of this RMT 
model, there actually exists a forth regime (Anderson 
strong localization regime) that in the present context 
does not have a semiclassical analog. 
As explained in M the generic hierarchy is 5x N u <C 
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